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Abstract 

I present a unified calculation of soft-gluon corrections to hard-scattering cross sections 
through next-to-next-to-next-to-leading order (NNNLO). Master formulas are derived, 
from a threshold resummation formalism, that can be applied to total and differential 
cross sections for hard-scattering processes in hadron colliders. I also present numerical 
results for charged Higgs production at the LHC where these corrections are large, and 
for top quark production at the Tevatron where these corrections greatly reduce the scale 
dependence of the cross section. 



1 Introduction 



Calculations of hard-scattering cross sections become very complicated as one moves from the 
lowest order to higher-order corrections. Current theoretical approaches include a variety of re- 
summations and fixed-order calculations, some through next-to-next-to-leading order (NNLO) 
PP. These tools can greatly improve next-to-leading order (NLO) calculations because higher- 
order corrections reduce the scale dependence and increase theoretical accuracy. 

The QCD corrections can be separated into hard, soft, and virtual parts, corresponding to 
contributions from energetic, soft, and virtual gluons, respectively. The soft-gluon corrections 
are an important component of the total result and, in some schemes and kinematical regions, 
e.g. threshold, they are numerically dominant. There is a universality in the form of these soft- 
gluon corrections, as can be clearly seen from the techniques of threshold resummations, which 
formally resum the soft-gluon contributions to all orders in the strong coupling j2j El IH 13 El 13 EJ ■ 
Although resummed calculations are prescript ion- dependent (see discussion in Ref. finite- 
order expansions of resummed cross sections are not, and they have provided us with many cross 
sections with NNLO soft-gluon corrections (for a review see Ref. PU|). Examples include W- 
boson production at large transverse momentum [TT], direct photon production |12j . top quark 
production 01131111, bottom and charm quark production [T3], heavy quark electroproduction 
[To] and photoproduction |T7j, charged Higgs production [TH], and jet production [TUj . 

The NLO soft-gluon corrections are typically a very good approximation to the exact NLO 
corrections near threshold. The NNLO soft-gluon corrections can be numerically significant 
and they invariably improve the theoretical calculation by stabilizing the dependence of the 
cross section on the factorization and renormalization scales, which are arbitrary energy scales 
in the theory. Thus it is worthwhile to provide a unified approach for the calculation of these 
and even higher-order soft-gluon corrections to hard-scattering processes in hadron colliders. 
It is important to note that new particles, such as in supersymmetry [20] or Higgs physics 
|21j . or particle production via new interactions, such as top production via flavor- changing 
neutral currents [22] , will likely be discovered near threshold where the soft-gluon corrections are 
dominant, therefore soft-gluon calculations are relevant to more than just pure QCD processes. 

The calculation of hard-scattering cross sections in hadron-hadron or lepton-hadron colli- 
sions can be written as 



where a is the physical cross section, <f>f/hi is the distribution function for parton / carrying 
momentum fraction Xi of hadron h{, at a factorization scale /xp, and hr is the renormalization 
scale. The parton-level cross section is denoted by &, and s, ti are standard kinematical invari- 
ants formed from the 4-momenta of the particles in the hard scattering. In a lepton-hadron 
collision we have one parton distribution (i — 1) while in a hadron-hadron collision we have 
two, % = 1,2. The partonic processes are of the form 



where / [I] represents a parton [lepton], F represents an observed system in the final state, 
and X any additional allowed final-state particles. For example, F can represent a pair of 
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heavy quarks, a single heavy quark, a jet, a photon, a Higgs boson, a pair of squarks, etc. 
Then s — (pi + P2) 2 - In single-particle-inclusive (1PI) kinematics we identify one particle 
F with momentum p, and define the kinematical invariants t\ = {pi — p) 2 , t 2 = (p% — p) 2 
(also commonly denoted by t and u, respectively). In pair-invariant-mass (PIM) kinematics we 
identify a pair of particles (such as a heavy quark antiquark pair) with invariant mass squared 
Q 2 . We note here that a and a are not restricted to be total cross sections; they can represent 
any differential cross section of interest. 

In general, a includes plus distributions T>i(x t h) with respect to a kinematical variable x t h 
that measures distance from threshold, with / < 2n — 1 at nth order in a s beyond the leading 
order. These are the soft corrections. The virtual corrections multiply delta functions S(xth)- 
In 1PI kinematics, x t h is usually denoted as s 4 (or s 2 ), s 4 = s + ti + t 2 — J2 m2 (the sum is 
over the masses squared of all particles in the process), and it vanishes at threshold. The plus 
distributions are of the form 
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where M 2 is a hard scale relevant to the process at hand, for example the mass m of a heavy 
quark, the transverse momentum p? of a jet, etc. The distributions are defined through their 
integral with any smooth function, such as parton densities, by 
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In PIM kinematics, with Q 2 the invariant mass squared of the produced pair, xth is usually 
called 1 — x or 1 — z, with z = Q 2 /s — > 1 at threshold. Then the plus distributions are of the 
form 
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defined through their integral with any smooth function by 



dzf(z) 
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The highest powers of these distributions in the nth-order corrections are the leading logarithms 
(LL) with / = 2n — 1, the second highest are the next-to-leading logarithms (NLL) with I = 
In — 2, the third highest are the next-to-next-to- leading logarithms (NNLL) with I = In — 3, 
the fourth highest are the next-to-next-to-next-to-leading logarithms (NNNLL) with I = In — 4, 
etc. These logarithms can be in principle resummed to all orders in perturbation theory. By 
now there are several processes for which NLL resummations and NNLO-NNLL results (i.e. 
the NNLL terms 1 at NNLO) have been presented [TTIll2"3] . 

lr rhe counting of logarithms is different in the exponent and in the fixed-order expansions; for example a 
term that is NNLL in the fixed-order expansion, as described here, may be NLL in the exponent |13j . 
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In this paper I present master formulas for the NLO, NNLO, and NNNLO soft-gluon cor- 
rections for processes in hadron-hadron or hadron-lepton collisions. These processes can be of 
QCD, electroweak, Higgs, or supersymmetric origin at lowest order. Results on the NLO and 
NNLO corrections have been presented before (22! but the notation here is somewhat different 
to facilitate the calculation of the NNNLO corrections. The NNNLO results and their appli- 
cations to top quark and charged Higgs production are new. In the next section, I present 
a threshold resummation formula from which high-order expansions are derived. In Sections 
3 and 4 are presented master formulas for the NLO and NNLO soft corrections, respectively, 
that arise from the expansion of the resummation formula. The formulas are given in the MS 
scheme (see (22] for results through NNLO in the DIS scheme), and cover both 1PI and PIM 
kinematics. In Section 5, I present a master formula for the NNNLO corrections. In Sections 
6 and 7 applications to charged Higgs production at the LHC and top quark pair production 
at the Tevatron, respectively, are discussed. Conclusions are given in Section 8. Some long 
expressions for terms in the NNNLO master formula are collected in an Appendix. 



2 Soft-gluon corrections from threshold resummation 

We begin with a brief review of the threshold resummation formalism. Threshold resumma- 
tion follows from factorization theorems for hard-scattering cross sections. One can write the 
hadronic cross section as a convolution of parton densities with a parton-level cross section. 
This can be further refactorized into functions associated with soft and collinear gluon emis- 
sion from the incoming partons and any outgoing partons or jets, a function associated with 
noncollinear soft gluon emission that involves the color structure of the hard scattering, and 
a short- distance hard-scattering function. The renormalization group properties of these func- 
tions result in the exponentiation of the soft-gluon contributions thus providing the resummed 
cross section [7] . For a review see Ref . . 

The resummation of threshold logarithms is carried out in moment space. We define mo- 
ments of the partonic cross section by &(N) = J dz z N ~ 1 a(z) (PIM) orbya(iV) = J^ds^/s) e~ Ns4 ^ s a(s4) 
(1PI), with N the moment variable. The logarithms of iV exponentiate. The resummed partonic 
cross section in moment space is then given by 
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The sums over % run over incoming partons: in hadron-hadron colisions we have two partons 
in the initial state, so % = 1, 2; in lepton-hadron collisions we have one parton. The sum over j 
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is relevant if we have massless partons in the final state at lowest order. We note that we have 
suppressed all gauge-dependent terms because these terms cancel out explicitly. 

Equation (j2.1j) is actually valid for both 1PI and PIM kinematics with appropriate definitions 
for Ni and Nj. In 1PI kinematics N{ = N(—ti/M 2 ) for incoming partons i, and Nj = N(s/M 2 ) 
for outgoing partons j; here M 2 is any chosen hard scale relevant to the process at hand. In 
PIM kinematics iVj = Nj = N. Also note that N = A^e 7B , with 7^ the Euler constant. The 
various exponents above are known at most to three loops. Below we give explicitly only the 
one-loop (and some two-loop) expressions that will be needed in our applications to charged 
Higgs production at NLL accuracy and top quark production at NNLL accuracy. Some two-loop 
and three- loop results can be found explicitly in and [SHIES]- We note however that the 
two-loop and higher-loop results for pro cess- dependent functions, such as the soft anomalous 
dimensions r$, have to be calculated explicitly for a specified partonic process. 

The first exponent in Eq. (|2.1|) is given in the MS scheme by 

EHNi) = - / dz I ^ -yAi (a s (Xs)) + v t [a s ((l - z) 2 s)} , (2.2) 

Jo I — z {J(i-z) 2 A L J J 

with Ai(a s ) = A^as/n + Af\a s /n) 2 + Af\a s /n) 3 + ■■■. Here A^ = Q with C t = C F = 
(N 2 — 1)/ (2iV c ) for a quark or antiquark and Cj = Ca = N c for a gluon, with iV c the number of 
colors, while A^ = CiK/2 with K = Ca (67/18 — 7r 2 /6) — 5n//9 [2Zj, where n/ is the number 
of quark flavors. Also Ui = (a s /7r)^- + (a s /7r) 2 v^ + (ag/ir) 3 ^ + • • •, with = C*. 
The second exponent is given by 



E' fj (Nj) = ( dz Z -^- / * ^Aj (a s (Xs)) - Bj [a s ((l - z)s)\ - v 3 k((l - zft 

JO 1 — Z I J (I— z) 2 A L 

(2.3) 

where Bj = (a s /n)Bf + (a s /Tr) 2 Bf + {a s /ixfBf ) + ■ • • with 5« = 3C F /A and = fa/ A, 
where (3q is the lowest-order /3-function, @ = (IICa — 2n/)/3. 

In the third exponent 7^ is the moment-space anomalous dimension of the MS density 

4>i/i. Note that the iV-independent part of the one-loop 7^ is the same as 7^, the one-loop 
parton anomalous dimension, given by 7^ = 3CV/4 and 7^ = /3q/4 for quarks and gluons, 
respectively [T3] . 

The (3 function in the fourth exponent is given in the Appendix. The constant d as =0,1,2 
if the Born cross section is of order a°, a], a 2 , respectively. 

fffifj are the hard-scattering functions for the scattering of partons fi and fj, while S' 1 *' 
are the soft functions describing noncollinear soft gluon emission. We use the expansions 
H = a^H® + (a d s «° +1 /ir)HM + ( a d s ^+ 2 /tt 2 )H^ + (af-+ 3 /n 3 )H^ + ■ ■ ■ and S = S<® + 
(a s j V) + (a s j n) 2 + (a s / n) 3 + • ■ -. Note that both H and S are matrices in color space 
and the trace is taken. At lowest order, the trace of the product of the hard matrices H and soft 
matrices S reproduces the Born cross section for each partonic process, a B = af as tr[if(°) S^]. 
The evolution of the soft function follows from its renormalization group properties and is given 
in terms of the soft anomalous dimension matrix Ts 00 121] • ^ n processes with simple color 
flow is a trivial lxl matrix while in processes with complex color flow an appropriate choice 
of color basis has to be made. For quark- (anti) quark scattering, is a 2 x 2 matrix 0I2E]; for 
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quark-gluon scattering it is a 3 x 3 matrix j^J; for gluon-gluon scattering it is an 8 x 8 matrix 
0. For the discussion below we expand T s as T s = (a s /7r)rg + (a; s /7r) 2 r<?' + (ag/ny^T^ + ■ ■ ■. 
The pro cess- dependent soft anomalous dimension matrices have by now been presented at one 
loop for all 2 — > 2 partonic processes; a compilation of results is given in |24j . They can 
be explicitly calculated for any process through the calculation of eikonal vertex corrections 
using the techniques and results in Refs. [5J El 121] • Some work has been done on two-loop 
calculations of these anomalous dimensions [2H], and furthermore the universal components of 
these anomalous dimensions for quark-antiquark and gluon-gluon initiated processes have been 
extracted from the NNLO results for Drell-Yan and Higgs production as detailed in Ref. [2*3*] . 

The exponentials in the resummed cross section can be expanded to any fixed order in a s and 
then inverted to momentum space to provide explicit results for the higher-order corrections. 
A fixed-order expansion avoids the problems with infrared singularities in the exponents and 
thus no prescription is needed to deal with these in our approach (see discussion in Ref. 0). 



3 NLO master formula for soft-gluon corrections 

We first expand the resummed formula in Eq. (j2.1j) to next-to-leading order and present a 
master formula for the NLO soft-gluon corrections in the MS scheme and 1PI kinematics: 
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where a B is the Born term, 
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with C q = Cp and C g = Ca, and c 2 is defined by c 2 = c% + T 2 , with 
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denoting the terms involving logarithms of the scale, and 
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denoting the scale-independent terms. We remind the reader that the sums over i run over 
incoming partons and the sums over j run over any massless partons in the final state. Note 
that not all the NLO corrections are proportional to the Born term; only the leading logarithms 
and terms involving the scale are. The function A c is process-dependent and depends on the 
color structure of the hard-scattering. It is defined by 

A c = tr (H^T^S® + HV>SV)r$>) . (3.5) 

With regard to the 5(xth) terms, we split them into a term c\, that is proportional to the 
Born cross section, and a term T{ that is not. C\ = d{ + T\, with 
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denoting the terms involving logarithms of the scale. T\ and Tf do not involve the factorization 
and renormalization scales. Note that T% and Tf are virtual terms and cannot be derived from 
the resummation formalism, but they can be read off by matching to a full NLO calculation 
for any specified process. 

In PIM kinematics we simply replace by 1 — z, set s = M 2 , and delete all ln(— ti/M 2 ) 
terms from the above expressions. The same should be done for the NNLO and NNNLO results 
that follow. 

As shown in Ref. j22] the NLO master formula passes a number of tests. Its predictions 
agree with NLO soft-gluon results for all processes where those results are already available. 
Also, the renormalization and factorization scale dependence in the physical cross section (after 
convoluting the partonic cross section with the parton distributions) cancels out explicitly, i.e. 
da/dfip = and da/du jR = at NLO. 



4 NNLO master formula for soft-gluon corrections 



At next-to-next-to- leading order, the expansion of Eq. (|2.1|) . with matching to the NLO soft- 
plus- virtual result, Eq. (j3.1|) . gives the NNLO soft corrections in the MS scheme and 1PI 
kinematics: 
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a 2 ^ 2 
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We have used the definitions 
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In PIM kinematics simply replace S4 by 1 — z, set s = M 2 , and delete all ln(— ti/M 2 ) terms. 

The quantities d aa , /3o, and K have all been defined in Section 2, and (2 = k 2 /6, (3 = 
1.2020569 • • •. Also C3, C2, ci, Tf, and A c have been defined in Section 3. The virtual terms 
and cannot be derived from resummation. A separate calculation is needed for each process 
to derive those. However, all the scale-dependent terms in R^ (R^ is scale-independent) can 
be derived and are given explicitly in Eq. (12) of Ref. |10j . 

As shown in Ref. j2H] , the NNLO master formula passes many rigorous tests. It reproduces 
the NNLO soft-gluon results for all processes where these results are known. Also at NNLO the 
renormalization and factorization scale dependence in the physical cross section cancels out. 

The master formula can in principle provide all the soft corrections at NNLO for any process. 
In practice, the accuracy which we can attain depends on whether the one-loop is known 
(in which case we can attain NLL accuracy; I^ is known for all 2 — > 2 processes); whether 

furthermore the NLO virtual terms are known (NNLL accuracy); and whether the two-loop 

(2) 

T s is known (NNNLL accuracy). Most current results are known to NLL or NNLL accuracy. 

(2) 

Note that T s is only known for the simplest cases of Drell-Yan and Higgs production where 

the color structure is trivial |2Sj- However, it was shown in [Hj that the contributions of T s 

(2) 

can be small so that effectively NNNLL calculations can be made in some cases even when T s 
is not fully known. 

5 NNNLO master formula for soft-gluon corrections 

At next-to-next-to- next-to-leading order, the expansion of Eq. (|2.1|) . with matching to the NLO 
and NNLO soft-plus- virtual results, gives the NNNLO soft-gluon corrections in the MS scheme 
and 1PI kinematics 
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The quantities 0^, /3 , an d have been defined in Section 2, while c 3 , c 2 , Ci, Tf, and A c have 
been defined in Section 3. Also C^, C^j, F c , G c , (2, and £ 3 have been defined in Section 4, 

and C 4 = vr 4 /90. The expressions for ft, X 3 , X 2 , X 1; X , Y; (1) , K^, f/ 0) , y/ 0) , M c , iff, if 2 c , 
and K3 are given in the Appendix. 

Note again that in PIM kinematics we simply replace S4 by 1 — z, set s = M 2 , and drop the 
terms with ln(— ti/M 2 ) in the above formula and in all the expressions given in the Appendix. 

This NNNLO master equation gives the structure of the NNNLO soft corrections and can 
provide the full soft corrections explicitly if all the two-loop and three-loop quantities are 
known. Therefore for processes with non-trivial color structure we are currently limited to 
NLL or NNLL accuracy, as the applications to charged Higgs and top quark production in the 
next two sections illustrate. The structure of the corrections as presented here can be useful 
for checking future calculations if and when such three-loop quantities become available. Also 
note that scale logarithms and Q constants will be kept as appropriate at subleading logs as 
explained in the next two sections. 

When the color structure of the hard scattering is simple, i.e. when H, S, and T$ are simply 
lxl matrices, then the above expressions can be simplified. We can then easily absorb A c into 
c 2 in Eq. ()3.1|) . and Tf into c±, by redefining c 2 and c\. Then all the terms are proportional to 
<j B in Eqs. (j3.1|) . (j4.1j) . and ()5.1|) . We will see this explicitly in the two applications in the next 
two sections. 
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6 Charged Higgs production via bg — > tH 

Charged Higgs production is a process of great interest at the LHC. The charged Higgs boson, 
if discovered, would be an unmistakable sign of new physics beyond the Standard Model [2*T] . 
A promising channel of discovery is associated production with a top quark via bottom-gluon 
fusion for which SUSY and QCD radiative corrections have been calculated [IJUJ EU EH!- 
NLO and NNLO soft-gluon corrections to this process were recently studied in ^H] where the 
corrections were found to be large, especially for a very massive charged Higgs. 

We now apply our NNNLO master formula to charged Higgs production in the Minimal 
Supersymmetric Standard Model (MSSM) via bottom gluon fusion, a process with simple 
color flow, at NLL accuracy We study the process b{p h ) + g(p g ) — > t(p t ) + H~{p H -) in 1PI 
kinematics and define the kinematical invariants s = (p& + p g ) 2 , t = (pb — p t ) 2 , u = (p g — p t ) 2 , 
and S4 = s + t + u — m 2 — m 2 H - , where m t is the top quark mass and m H - is the charged Higgs 
mass (we take the bottom quark = in the kinematics ^K|) • 

The NLO coefficients of Section 3 here take the values C3 = 1{Cp + Ca) and c 2 = c 2 + T 2 
with c£ = -{C F + C A ) ln(/4,/m|_) and 

T 2 = 2Relf -C F -C A - 20 In ( ~ U + ) - 2C A In ( Zl±<A 



m H- ) \ m H- 



' 1 

m 



- (C F + C A ) In ( — J-J . (6- 1 ) 

where Rer^ denotes the real part of the one- loop soft anomalous dimension [18] . 
r (i) n , f-t + m 2 \ C A (-u + m 2 \ C A 

Note that, due to the simple color structure of this process, Ts is simply a 1 x 1 matrix. Here, 
as described in the last paragraph of the previous section, we have aborbed the term 2Rer^ 1 ' ) , 
which arises from A c , into T 2 . Also, 



—u + m 2 . 



Oln "V g " +Ca1h 





In 



' /4 



(6.3) 

The various terms from the Appendix used in the NNNLO corrections here take the values 
X 3 = /W12, X% = ((3 /8)c 3 \n(ii 2 R /m 2 H _), 

^ -Mi) - (c ^Mt ) • m 

X} = j (2 c 3 , (6.5) 

2 

where X% denotes the scale logarithm terms in X 2 , Xf denotes terms involving squares of the 
scale logarithms in Xi, and X^ denotes the Q terms in X%. 
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Figure 1: The total cross section for charged Higgs production at the LHC. 



The NNNLO-NLL corrections are then given by 



a 



(3) 



a 



7T" 



+ 
+ 

+ 



c 3 (c£) 2 + 2 c 3 T 2 + £ c z 3 d{ - ( 2 - 4 c£ X 3 + 2 c 3 X£ 



P 5 (s 4 ) + 



,2 u 



5 2 5 

- c 3 c 2 - - c 3 X 3 



VAsA 



3 1 

- c 3 c 2 q + - [c 2 



3 ( 2 c 2 c 2 + - C3 c 3 



+ y c 2 c 3 x 3 + 3 x£ - 2 C3 ( xf 2 + xf; 



VAsA 



{4Y c? -tictdi-- C 2 c 3 ((c^ + 2 T 2 Ai) + 5 Cs c 2 

+ 12(2^X3-5(203X^-2^ (xf'+xf)] Px^)} . 



(6.6) 



Because we absorbed A c into c 2 the corrections take a simple form, simply multiplying the 
Born term a B . Note that consistent with a NLL calculation we include all T>^ (LL) and T>^ 
(NLL) terms. In addition, we calculate all scale logarithms at NLL accuracy. This means that 
for coefficients of ln*(/x 2 /m 2 ) we include the most singular plus distribution and the next-most- 
singular one ^31- Thus, we also include all scale logarithms in the D 3 terms, the cubed and 
squared scale logarithms in the T> 2 terms, and the cubed scale logarithms in the T>\ terms. 
With respect to the subleading Q terms that arise from inversion from moment to momentum 
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Figure 2: The .fT-factors for charged Higgs production at the LHC. 



space, we include only those that we can calculate exactly (for a discussion of the numerical 
effects of such terms see Ref. [9 ). Thus we include all Q terms in the V 3 and V 2 terms, and 
all Q terms multiplying scale logarithms in the T>\ term. 

In Fig. plot the cross section versus charged Higgs mass for pp collisions at the LHC 
with v^S* = 14 TeV using the MRST2002 approximate NNLO parton distribution functions 
|3"4"j with the respective three-loop evaluation of a s . We set the factorization scale equal to the 
renormalization scale and denote this common scale by fi. We show results for the LO, NLO- 
NLL, NNLO-NLL, and NNNLO-NLL cross sections, all with a choice of scale fi = m^- . In 
the calculation we choose a value tan f3 = 30; here tan (3 is the ratio of the vacuum expectation 
values of the two Higgs doublets in the MSSM. It is straightforward to calculate results for 
any other value of tan (3, since the only dependence on (3 is in a factor m 2 tan 2 (3 + m 2 cot 2 (3 
appearing in the Born term. The cross sections span over two orders of magnitude in the 
mass range shown, 200 GeV < m H - < 1000 GeV. The NLO, NNLO, and NNNLO threshold 
corrections are positive and provide a significant enhancement to the lowest-order result. We 
note that the cross sections for the related process bg — > iH + in the MSSM are exactly the 
same. 

In Fig. |^1 we plot the .fT-factors, i.e. ratios of cross sections at higher orders to the LO 
result, to better show the relative size of the corrections. The NLO-NLL / LO curve shows 
that the NLO soft-gluon corrections enhance the LO cross section by approximately 25% to 
50% depending on the mass of the charged Higgs. As expected the corrections increase for 
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higher charged Higgs masses as we get closer to threshold. The NNLO-NLL / LO curve shows 
that if we include the NNLO threshold corrections we get an enhancement over the LO result 
of approximately 35% to 70% in the range of masses shown. Again the enhancement increases 
with charged Higgs mass, as expected. Finally, the NNNLO-NLL / LO curve shows the further 
enhancement that the NNNLO soft-gluon corrections provide, approximately 45% to 80% over 
LO. We note that the NNNLO soft-gluon corrections are quite significant in this case. This is 
not necessarily typical for other processes, but happens here because of the very massive final 
state. Another process where NNNLO soft-gluon corrections are known to be big is inclusive 
hadron production at high transverse momentum |35j . 



7 Top quark production via qq 



tt 



The study of the top quark is important in understanding the electroweak sector and searching 
for new physics. The top quark is now being actively studied at Run II at the Tevatron 
[Hot l3*7| I3*H] . Theoretically, the production cross section has been studied to NNLO and NNNLL 
in both 1PI and PIM kinematics 14j. The corrections are moderate and they substantially 
decrease the scale dependence of the cross section. Transverse momentum distributions are 
also known to NNLO-NNNLL ^3], while rapidity distributions have been presented in Ref. 



Here we apply our NNNLO master formula to top quark pair production via quark-antiquark 
annihilation, which is the dominant partonic subprocess at the Tevatron. We study the channel 
q(Pa) + q(pb) t(pi) + t(p2) in 1PI kinematics, and define the kinematical invariants s = 



(Pa+Pb) 2 , h = (Pb-PIY 



m 



t ■ 



Ul = (p a - Plf 



m 



t ) 



and S4 = s + t% + u\, with m t the top 



quark mass. For this process, which has complex color flow, H, S, and T$ are 2x2 matrices. 
However, has a particularly simple form in a singlet-octet color basis (the only non-zero 
element is jUlEI]) so that simplifications arise. 

The NLO coefficients of Section 3 here take the values c 3 = 4C F and c 2 = c 2 + T 2 , with 
-2C F ln(fi 2 F /m 2 ) and 



u 
("2 



T, 



2 Rer^ 22 - 2C F - 2C F In 



t-iU 



iui 



ml 



- 2C F In 



(7.1) 



where Rer^ 22 denotes the real part of the 22 element of the soft anomalous dimension matrix 
0i, 



5,22 



C7 



41n(^ 
ti 



+ 9a 

2 



-31n(^ 

U 



In 



Uu Xi 



+ Lp + in 



(7.2) 



with L/3 = [(1 - 2m t 2 /s)//3][ln((l - (3)/(l + (3)) + in] and (3 = ^/1- 4m 2 t /s. Note that even 
though we have nontrivial color matrices for this process, due to the simple form of the H 

To be precise, 



matrix here we have aborbed the term 2Rer^ 22 , which arises from A c , into T; 



A c = ((r B /a 2 ) 2 ReT { s\ 2 . Also, c x = d[ + T\ with 



In 



In (^)+^ln 

\m 2 t 2 \m? 



(7.3) 
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qq -> tt at the Tevatron S =1.96 TeV m=175 GeV 
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Figure 3: The scale dependence of top production in the qq channel at the Tevatron. 



and T\ as given by the 5{s^) terms in Eq. (4.7) of Ref. [lOj (we have effectively absorbed 
Tf into ci; for details see [II]). Also, F c = (a B /a 2 s ) [4(R e rg 2 ) 2 + 4rg ) 12 rg ) 21 ] with r£l 2 = 



(Cf/Ca) ln(wi/ti) and r^ 21 = 21n(-ui/ti) the 12 and 21 elements of the soft anomalous dimen- 
sion matrix. 

The various terms from the Appendix used in the NNNLO corrections here take the values 
X 3 = fa ca/12, X 2 = -(A)/4)T 2 + (A)/8)c 3 H^t/ml) + c 3 K/i, X£ = (ft/ 8)c 3 ln^/mf), 



Xf 



Xf 



4 

/3o 



mi 



4 



mi 



T T 2 ln (&)+C F Kbx(!£ 
4 \mi \mf 



(7.4) 

(7.5) 
(7.6) 



+ 



0o 



-C F - C F In 



16 

'hui 



mi 



mi 



In 



2 / /fl 



(7.7) 



15 



Ik 



C3C3 - ^(2^2 + y C2C3 In 



8 



Y * ~ ~ Cf 48 



/4 



,/4 



(7.8) 
(7.9) 
(7.10) 



. .71 /- 

where Xf denotes the terms with scale logarithms to the n-th power in Xi, Xq denotes the Q 
terms in Xo, Yj 1 '^ ^ denotes terms involving squares of the scale logarithms in Y^\ and Y^ 0,At ) 
denotes terms involving cubes of the scale logarithms in Yj°\ 
The NNNLO-NNLL corrections are then given by 



a 



(3) 



a 



R) 



TV' 



-c 3 3 V 5 ( Si ) + 



5 2 5 

— Co Co — — Cq 

8 3 2 



+ 
+ 

+ 



c 3 c\ + x - c 2 ci - G cjj + (/So - 4 c 2 )X 3 + 2 c 3 X 2 + 4 c 3 rg } 12 rg 21 



^ 3 (s 4 ) 



|c3C^ Cl + |c 3 T 2 c^ + ^ (c\- T|) -3C 2 c2c 2 + |c3C 3 -3c^X3 + ^C2C3X 3 



(3 T 2 - ft) X£ + 3 c^X 2 - - c 3 (xf + xf 2 + x{) + 6 c^ rg{ 2 r« 



^ 2 (s 4 ) 



c 3 (cO 2 + (c^) 2 Cl + 2 T 2 c^ c? - C2 c 2 ci - (2 c 3 c 2 + 5 Cs c 2 c 2 + ^ Cl c 3 

-YC2C3 + (-20C3C 3 + 12C 2 c 2 )X 3 + 2c^X 2 ^-5C 2 c 3 X 2 + (/?o-2c 2 ) (xf+X$) 
-2<$X? + c,(xf +X C C 
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+ 2y^ 2 )-ioc 2C3 rgi 2 rg 21 



£>i(s 4 ) 



+ 



2^«) 2 



— (4 c 3 c^ + Cs c 3 c^ + 2 Cs c 3 (c 2 - T|) + - C| c 3 c^ - (2 c 3 c 2 c£ 

C 2 c 3 c, 1 T\ - ~ C 2 (c 3 - T|) + ( -8 C 3 4 + 3 C2 cl) X 3 + (4 Cs c 3 - 3 ( 2 T 2 ) X£ 

- (Xf 2 + Xf) + c, 1 (x£ 2 + X C ) + 2 Ff^ 3 ) 



-3C 2 ^X 2 + C 2 c 3 (xf + xf 

-ec^rg^rgy d q (s 4 )} - 



(7.11) 



Because we have absorbed v4 c into c 2 , and T x c into ci, the corrections take a simple form, 
simply multiplying the LO term a B . Note that consistent with a NNLL calculation we include 
all T>5, T>4, and £> 3 terms. In addition, we calculate all scale logarithms at NNLL accuracy. This 
means that for coefficients of In 1 (/i 2 /m 2 ) we include the most singular plus distribution and the 
two next-most-singular ones ^3]. So we also include all scale logarithms in the T> 2 terms, the 
cubed and squared scale logarithms in the T>\ terms, and the cubed scale logarithms in the 
T>q terms. With respect to the subleading Q terms that arise from inversion from moment to 
momentum space, we include only those that we can calculate exactly (see Ref. 0). Thus we 
include all Q terms in the T> 2 and T>\ terms, and all Q terms multiplying scale logarithms in 
the T>q term. 
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In Fig. |3]we plot the scale dependence of the cross section for qq —>■ tt at the Tevatron Run II 
for a top quark mass rrtt = 175 GeV. We set the factorization scale equal to the renormalization 
scale and denote this common scale by \i. Again we use the MRST2002 approximate NNLO 
parton distribution functions JS] (results using the CTEQ6M distributions |4T] are similar 
|14j). We plot a large range in scale and see that the higher-order soft-gluon corrections greatly 
decrease the scale dependence of the cross section. The NNLO-NNNLL and the NNNLO- 
NNLL curves are relatively flat. The NNNLO-NNLL result displays very little variation and it 
approaches the scale independence expected of a physical cross section. 



8 Conclusions 

In this paper, I presented a unified approach to calculating the NNNLO soft-gluon corrections 
for hard-scattering processes in hadron-hadron and lepton-hadron collisions in the MS scheme in 
either 1PI or PIM kinematics. The master formulas given in the paper allow explicit calculations 
for any process, with either simple or complex color flows, keeping in general the factorization 
and renormalization scales separate and the color factors explicit. Detailed results, illustrating 
the use of the master formulas, were given to NLL accuracy for charged Higgs production 
via bottom-gluon fusion at the LHC, and to NNLL accuracy for top quark production via 
quark-antiquark annihilation at the Tevatron. 

The NNLO and NNNLO corrections increase theoretical accuracy and diminish the de- 
pendence on the factorization and renormalization scales, and thus are essential in further 
calculations of QCD corrections for many hard-scattering processes, both in their own right 
and as backgrounds which may be particularly important in searching for the Higgs boson 
and supersymmetric particles, as well as other processes that signal new physics beyond the 
Standard Model. 



Appendix 

The (5 function is given by 

P(a.) = tiding /dyL = -/W(^) - Aa^/(4vr) 2 - P 2 a 3 J(^) 3 + ■ • • , (A.l) 
where g 2 = 4vra s , with fa = {UC A - 2n/)/3, ft = UC\/?> - 2n f (C F + 5C A /3), and 

ft = ^ + (d - ^C F C A - ^d) n, + (§C F + IC A ) n) . (A,) 

Note that 
a s {n) = a a (jj, R ) 1 - 

47T 

+ •••]. (A.3) 

The various quantities used in the NNNLO expression, Eq. (|5.1|) . are given by the following 
expressions. For the quantities X 3 , X 2 , X\, X , we have 

* = (A.4) 

3 
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a s (fi R ) In (4) + ^«^W In 2 (4) - A«20*h) In (4) 
\PrJ 167r \PrJ 167r \»R/ 



V " 7 '°r It, / ^ _ n V-^0 R (1) 



K 



(A.5) 



Xi = c 2 ci - C2 c 3 c 2 + Cs C3 + ^ C2 c 3 - E C i^ 2 ~ C d] > 



(A.6) 



and 



X = ^ Cl lnf41+ — 
4 1 M 2 16 



+ f (sea - |c 2 r 2 + fc^ m (§) + E af* - E y$ m (§) 



A) 



7< W " C7* In 



+ E^4 ln 



2 "* \M 2 J " I M 2 



The quantity i?'^ at the end of the above equation stands for the virtual two-loop corrections 
in Eq. (@~U) minus £ terms and scale terms. To be precise 



d(2) d'(2) 1 1 2 C2 2 , 1^2 2 , /■ % 2 8q f flR 
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and 

+ 3H (o } (r^t) V)r« + 3ff<°>rg> W> (rg 
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